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Abstract 

We study a lattice QCD mixed action with overlap valence quarks on two flavours 
of Wilson maximally twisted mass sea quarks. Employing three different matching 
conditions to relate both actions to each other, we investigate the continuum limit 
by using three values of the lattice spacing ranging from a ~ 0.05 fm to 0.08 fm. A 
particular emphasis is put on the effect on physical observables of the topological 
zero modes appearing in the valence overlap operator. We estimate the region of 
parameter space where the contribution from these zero modes is sufficiently small 
such that their effects can be safely controlled and a restoration of unitarity of the 
mixed action in the continuum limit is reached. 
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1 Introduction 



The discretization of the QCD action on a lattice provides a framework for 
a quantitative description of the non-perturbative phenomena inherent to 
hadronic physics at low energies. The way of discretizing the continuum QCD 
action is, however, by no means unique and leaves a large number of choices. 
Although all lattice actions used in practise respect the principle of local gauge 
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invariance, they often break a number of continuum symmetries depending on 
the particular choice of the lattice action used. A most notable example is the 
explicit breaking of chiral symmetry on the lattice in order to avoid unwanted 
additional fermion excitations, the so-called doublers [TH3]- There exist, how- 
ever, lattice fermion actions which are constructed from Dirac operators that 
satisfy the Ginsparg- Wilson relation [1] and lead to an exactly preserved but 
lattice modified chiral symmetry |5]. The big disadvantage of these Ginsparg- 
Wilson type fermions is that they are computationally very expensive. 

A possible way to restore a desired symmetry -at least partially, since it will 
not by respected by the sea quark action - is the approach of so-called mixed 
actions, corresponding to the use of different discretizations of the Dirac op- 
erator in the valence and sea sectors. Of course, using such a mixed action it 
needs to be guaranteed that a well defined continuum limit, which corresponds 
to the universality class of continuum QCD, is reached eventually. This can be 
achieved by imposing an appropriate matching condition relating both actions 
to each other and which fixes the physical situation at each value of a non-zero 
lattice spacing. Hence, a well defined continuum limit can be reached. 

Although in this way mixed actions offer an advantage since the valence quarks 
can have improved symmetry properties with respect to the sea quark action, 
they also provide challenges. Even when imposing a matching condition, the 
sea quark and valence quark theories are not fully matched and differ by dis- 
cretization effects, which lead to unitarity violations in physical quantities. It 
is therefore an important question to understand the characteristics of a mixed 
action approach and to learn how and whether indeed the correct continuum 
limit is reached. In particular, it is an interesting and important question, 
whether and how the inherent mismatch of the eigenvalue spectra affect phys- 
ical observables and influence possibly the continuum limit. 

It is the goal of this work to address the above questions by studying a 
particular mixed action setup, employing maximally twisted mass (MTM) 
fermions [6l[7] in the sea and chirally invariant overlap fermions [8l[9] in the 
valence sector. As anticipated above, the mismatch of the spectra of the lattice 
Dirac operators employed in the sea and the valence actions will be a particu- 
lar concern. While the overlap Dirac operator exhibits exact chiral zero modes, 
for any non-zero value of the lattice spacing, the twisted mass Dirac operator 
does not, at least at lattice spacings used in this work. As discussed already in 
ref . [To] , the quantitative effects of these unmatched zero modes on physical 
observables will depend on the quark mass and the physical volume used in 
the numerical simulation. 

We will investigate the mesonic sector of the theory and study the contin- 
uum limit of light pseudoscalar meson observables by using three values of 
the lattice spacing ranging from a ~ 0.05 fm to 0.08 fm at a fixed lattice 
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size of L ^ 1.3 fm, and quark masses m^^(/i = 2 GeV) ~ 20MeV and 

m^^{fi = 2 GeV) ~ 40MeV. These values of the quark masses correspond to 
infinite volume pseudoscalar masses of approximately 300 MeV and 450 MeV, 
respectively. We also study finite size effects by considering three values of 
the lattice size between L ^ 1.3 fm and 2.0fm. The choice of this parameter 
region (quark mass, lattice size and lattice spacing) allows us to perform then 
a systematic study of our mixed action setup. 

Mixed actions are currently being used by various groups. In the case of valence 
overlap quarks, studies using domain wall sea fermions pT|[T2] or Wilson clover 
sea quarks [131 HI] have been carried out. Domain wall valence fermions on 
staggered sea quarks have also been extensively used [T5HT7]. Mixed actions 
in which variants of the sea action are used in the valence sector, such as 
Osterwalder-Seiler fermions on a twisted mass sea, have also been recently 
studied [T81 - I20] . The extension of chiral perturbation theory to the context 
of mixed actions [2T1423] provides a useful analytic control of the light quark 
mass dependence of several physical quantities in the regime of parameters 
(quark mass, lattice size and lattice spacing) where the effective theory can 
be applied. 

As said above, in this paper, we perform a systematic study of a mixed action 
with overlap valence fermions on a Wilson twisted mass sea, by focusing on 
light pseudoscalar meson observables. In Section 2 we introduce the mixed 
action setup by shortly reviewing the overlap and twisted mass formalisms 
and introducing the quark mass matching procedure, which is illustrated in 
the free theory. Section 3 describes the simulation setup and Section 4 presents 
the results of the continuum limit study of the mixed action for our light 
quark mass ensembles. Section 5 provides an interpretation and discussion of 
the results. In Section 6 we discuss the regime of heavier quark masses and 
finite size effects. Section 7 concludes and provides an estimate of the regime 
of parameters where the contribution in the valence sector of the topological 
zero modes on physical observables is small enough such that a proper control 
of the restoration of unitarity of the mixed action becomes possible. 



2 Mixed Action: Fermionic Actions, Matching and Illustration in 
the Free Theory 

2.1 Overlap Fermions 

Overlap fermions were introduced in 1997 by Neuberger [8l|9], who found a 
particularly simple form of a lattice Dirac operator that obeys the Ginsparg- 
Wilson relation [1]. This implies that overlap (ov) fermions preserve a lattice 
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version of chiral symmetry at any value of the lattice spacing [5] . The massless 
overlap Dirac operator is given byPI: 

DUO) = l{l-A{A^A)-'/'y (1) 

The kernel operator A reads: 

A = l + s-a£)wiison(0), (2) 

where s is a parameter which satisfies |s| < 1. The Wilson-Dirac operator is 
defined by: 

^Wiison(mo) = i (7/.(V; + V^) - aV;V^) + mo, (3) 

where mo is the bare Wilson quark mass and V^, V* are the forward and back- 
ward covariant derivatives, respectively. The massive overlap Dirac operator 
is given by: 

^ / CLin \ 
Do^{mo^) = f 1 £'ov(0) + mov, 

where mov is the bare overlap quark mass. Note that despite the appearance 
of the square root in eq. ([I]), the overlap operator is local in the sense of ex- 
ponentially decaying coupling strengths as a function of the distance between 
lattice points [25] . 

Operators preserving the Ginsparg- Wilson relation develop topological zero 
modes and fulfil the index theorem at finite lattice spacing [26]. Moreover, 
chiral symmetry allows to improve the theory at 0{a). Indeed, observables 
computed with overlap fermions are not affected by 0{a) lattice artefacts, 
provided that C(a)-improved interpolating operators are used (see e.g. [27]). 

The overlap quark mass rriov renormalizes multiplicatively with a renormali- 
sation factor Zm = where P stands for the pseudoscalar density. The 

decay constant /tt of a mass degenerate pseudoscalar meson, carrying a mass 
m^r, can be determined by means of the two-point pseudoscalar correlation 
function Cpp(t) via the following expression: 

/r = 7^l(0|P|vr)o.|. 
{mZ ) 



It is known that in practise, overlap fermions are demanding from the compu- 
tational point of view. With respect to Wilson type fermions the overall cost 
can grow by one to two orders of magnitude[l]- Moreover, the discontinuities 

^ For a review about overlap fermions, see e.g. f24J. 

^ We refer to |28] for a detailed comparison to Wilson twisted mass fermions. 
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of the overlap operator while changing topological sectors pose a problem to 
the standard Molecular Dynamics algorithm used in dynamical quark simula- 
tions. Modifications of the Molecular Dynamics evolution to tackle this prob- 
lem [2nH32] imply a further significant increase in the computational cost. A 
way to avoid the discontinuity problem is to modify the action in such a way 
that the topological charge remains fixed [33l|3l]. However, this procedure in- 
troduces additional 0{1/V) finite volume effects which need to be taken into 
account in the determination of physical observables. 



2.2 Wilson Twisted Mass Fermions 

Twisted mass (tm) fermions [6j are defined by adding a chirally rotated mass 
term to the Wilson-Dirac operator in eq. ^ as follows: 

-Dtm = -Dwilson("^o) + ^/^<?75T"3 , (4) 

where is the twisted mass parameter and is the third Pauli matrix acting 
in fiavour space. 

Wilson twisted mass fermions were introduced to address the problem of un- 
physically small eigenvalues of the Wilson-Dirac operator |6] . Furthermore, by 
tuning the mass rriQ to its critical value merit, a situation denoted by maximal 
twist, physical observables are automatically 0{a) improved [7] independently 
of the operator which is considered, implying that no additional, operator spe- 
cific improvement coefficients need to be computed. Detailed studies of the 
continuum- limit scaling in the quenched approximation [35] and with two dy- 
namical quarks [36j have demonstrated that, after an appropriate tuning pro- 
cedure to maximal twist, lattice artefacts indeed follow the expected 0{a?) 
scaling behaviour. 

A comprehensive simulation programme has been undertaken by the Euro- 
pean Twisted Mass (ETM) collaboration, including dynamical simulations 
with Nf = 2 [36H3H] and iVf = 2 + 1 + 1 [39] quark flavours. 

At maximal twist, the quark mass is given by the twisted mass /ig and, as 
in the case of overlap fermions, it renormalizes multiplicatively with a renor- 
malisation factor = 1/Zp. In contrast to standard Wilson fermions, an 
exact lattice Ward identity for maximally twisted mass fermions allows the 
extraction of the decay constant of the charged pseudoscalar meson from the 
relation 

which is free from the presence of renormalisation factors. 
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The twisted mass term in eq.(jl]) explicitly breaks parity and isospin symme- 
try, which are however restored in the continuum limit with a rate of O(a^) as 
shown in [7] and numerically confirmed in [351110] • The effect of isospin break- 
ing has been observed to affect in a substantial way only the neutral pion 
mass. Indeed, while the discretization effects in the charged pion are observed 
to be small, significant 0{a'^) corrections appear when studying the scaling 
to the continuum limit of the neutral pion [36]. Similar effects have not been 
observed in other quantities that are in principle sensitive to isospin break- 
ing but not trivially related to the neutral pion mass. These observations are 
supported by theoretical considerations detailed in [1T1, H2]. 

2.3 Matching Procedure 

In order to perform the continuum limit extrapolation, the physical situation 
(e.g. the quark mass and the lattice size) has to be kept fixed when changing 
the lattice spacing. In practise, this requires a matching procedure at each 
value of the lattice spacing. One possibility for this matching is to directly 
match the renormalized quark masses. Another option is to match a hadronic 
quantity that can be accurately determined with both regularisation and that 
depends in a significant way on the quark mass. In the light quark sector, the 
mass of the pseudo Nambu-Goldstone boson is a natural hadronic observable 
to perform this matching. In this work we follow this procedure, by defining 
the matching point in the following way: 

I mov = ''^!r™ I =^J.%'"^ ) ( 5 ) 

where the l.h.s is the mass of a pseudoscalar meson consisting of mass degen- 
erate overlap valence quarks of mass ruov = m^^^^ (the so called matching 
quark mass) in a background of iVf = 2 twisted mass sea quarks with mass 
fig = fj,^g^^. The mass m^™ in the r.h.s corresponds to the pseudoscalar meson 
made of twisted mass valence and sea quarks of mass fig = /x^®'^. 

To illustrate the matching procedure in the valence sector, we perform a con- 
tinuum limit scaling analysis of overlap and twisted mass fermions in the free 
theory. 

Illustration in the Free Theory 

We want to investigate the effects of the matching procedure of the quark 
mass in the free theory by using as a matching condition the pseudoscalar 
meson masses of maximally twisted mass fermions and of overlap fermions. 
The twisted mass is kept fixed to Nfig = 0.5, where is the number of lattice 
sites in the spatial directions. 
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Fig. 1. Tree-level examples : (a) the matching of twisted mass and overlap quark 
masses and (b) the mismatch between the twisted mass and overlap pseudoscalar 
decay constants at the matching point A^m^™ = Nm'^. 

We illustrate the matching procedure through an example with = 16. In 
this case, the twisted mass pseudoscalar meson mass is Nm}^ = 0.999959. 
Since the magnitude of 0{a'^) effects is in general different for different fermion 
discretizations, if we impose the condition of equal pseudoscalar meson masses 
Nm}^ = Nm'^, the bare quark masses will differ Nrriov 7^ Nfiq. This is shown 
in Fig. 1 (a) , where the dependence of the overlap pseudoscalar meson mass on 
the quark mass Nrriov is depicted. The value Nm^^^'^^ ~ 0.49994 leads to the 
same pseudoscalar meson mass as the value iV/Zg = 0.5 in the twisted mass 
case. 



Since the size of discretization effects depend on the particular observables 
under consideration, the matching of the pseudoscalar meson mass does not 
imply that other observables will also be matched. This is illustrated in Fig. 



1(b) which shows the difference of the pseudoscalar decay constants at the 
matching point, Nm}^ = Nm'^ . 

The difference A^/*™ — Nf^^ vanishes in the continuum limit at a rate of 
0(0^). We consider different lattice sizes, by changing in the range = 4 
to A^ = 64. By fixing the physical situation to Nfig = 0.5, the change in A^ 
introduces the scaling towards the continuum limit, which, in the free theory, 
corresponds to the limit A^ — 00 j^] For each lattice size A^, we determine the 
matching overlap quark mass Nm^^^^, such that the overlap and twisted mass 
pseudoscalar meson masses are equal. In Fig. [2] we illustrate the dependence 
of the difference Nf^ — A^/ps on 1/N'^. This difference indeed vanishes in 
the limit N ^ 00 and the leading discretization effects are of (9(1/A^^). A 
similar behaviour is observed for other observables such as the pseudoscalar 
correlation function at a fixed physical distance. 



^ More details about tree-level scaling tests of various fermion discretizations and 
their expressions for the quark propagators and correlation functions can be found 



m 
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0.00045 




Fig. 2. The difference, in the free theory, between the twisted mass and overlap 
pseudoscalar meson decay constants at the matching point Nni^ = Nm'^ , as a 
function of 1/N'^. The hmit corresponds to the continuum hmit. 

3 Simulation Setup 



In the following, we will investigate the just discussed example of free fermions 
in the case of including the interaction to the gluon fields. In table [1] we provide 
the list of dynamical simulations considered in this work. The ensembles were 
generated by the ETM collaboration [36] using the tree-level Symanzik im- 
proved gauge action and Nf = 2 flavours of maximally twisted mass fermions. 

The heavier ensembles labelled in table [1] with a subscript h, i.e. Bh, Ch 
and Dh, correspond to three values of the lattice spacing with a physical 
situation fixed to a lattice size L/tq 3 and a charged pseudoscalar meson 
mass m^™ro ~ 1.0, where tq is the Sommer scale |45j. This corresponds to 
a lattice size L ^ 1.3 fm and a non-perturbatively renormalized [IH] quark 
mass = 2 GeV) ~ 40MeV. In infinite volume this corresponds to a 

pseudoscalar meson mass m^r ~ 450 MeV. 

The lighter ensembles are labelled with a subscript i. In this case the pseu- 
doscalar meson mass is fixed to mTrTo ~ 0.8. This corresponds to a quark mass 
f^^^if^ — 2 GeV) ~ 20 MeV, which in infinite volume gives a pseudoscalar 
meson mass of around 300 MeV. The subscripts Vi^2.3 refer to ensembles with 
different spatial lattice extents. 

In the valence sector we use overlap fermions0- The gauge links entering in 
the covariant derivative of the kernel operator A in eq. are HYP smeared 
to reduce the condition number of A'^A. The mass parameter s is set to s = 
since this is observed to provide the best locality properties for the overlap 



More details about our overlap fermions setup are given in [i7] . 
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Ensemble 


/? 


{L/af X T/a 




^crit 




L/ro 


Bh 


3.90 


16^ X 32 


0.0074 


0.160856 


1.03 


3.0 


Be,Vi 




16^ X 32 


0.0040 




0.84 




Be,V2 




20^ X 40 






0.73 


3.8 


Be,V3 




243 X 48 






0.72 


4.6 


Ch 


4.05 


20^ X 48 


0.0060 


0.157010 


1.00 


3.0 


Ci 




20^ X 40 


0.0030 




0.83 




Dh 


4.20 


243 X 48 


0.0050 


0.154073 


1.04 


2.9 


Df 




24^ X 48 


0.0020 




0.82 





Table 1 

Summary of A'^f = 2 Wilson twisted mass ensembles considered in this work. We 
give the ensemble label, the values of the inverse coupling /3, the lattice volume 
X T, the twisted mass parameter a^g, the critical hopping parameter Kcrit = 
1/(8 + 2amcrit), the approximate values of the pseudoscalar meson mass and of 
the lattice size in units of the Sommer scale tq. The chirally extrapolated values 
of the Sommer scale are r^/a = 5.25(2) at /3 = 3.90, ro/a = 6.61(2) at (3 = 4.05 
and ro/a = 8.33(5) at (3 = 4.20, corresponding to lattice spacings a ^ 0.079 fm, 
a 0.063 fm and a 0.051 fm, respectively [36]. The subscripts h and £ in the 
ensemble name refer to the heavy and light quark masses, respectively, which are 
kept fixed when performing the continuum limit extrapolation. The subscripts 14^2,3 
label ensembles with different lattice sizes. 



operator 0. 

For each dynamical quark ensemble in table [H we produce all-to-all overlap 
quark propagators in a wide range of quark masses rriov, from the unitary (sea) 
light quark mass up to the physical strange quark mass. This provides a good 
handle on the matching procedure, to be discussed in the following section. 
Note that the lattice size L/r^ ^ 3 is outside the regime of applicability of 
chiral perturbation theory (or its generalisation to mixed actions). In this 
work, we instead rely on a direct study of the continuum limit of the mixed 
action in order to characterise the size of unitarity violations. 



^ Preliminary results were presented in [35] and we refer to a forthcoming publica- 
tion j39] for more details about this aspect. 
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4 Continuum Limit Scaling Tests and Different Matching Condi- 
tions 



In this section, we wiU discuss a number of continuum hmit scahng tests for 
the pseudoscalar decay constant. We wiU employ three different matching con- 
ditions to relate the overlap valence quarks to the maximally twisted mass sea 
quarks and the pseudoscalar decay constant will be extracted from two differ- 
ent correlation functions which differ by their sensitivity to chiral zero modes. 
We remark that we consider the matching condition employed to define the 
mixed action theory used. This is similar to improvement conditions where the 
resulting improvement coefficients define the improved actions and operators. 
In particular, there is then no uncertainty associated with a given match- 
ing condition. However, it is of course still important to investigate different 
matching conditions and study their effects, as we will do in this work. 

In the following we will use the shorthand notation "scaling test" for the 
continuum limit scaling test in the lattice spacing. As in the example of free 
fermions discussed above, we start with a naive scaling test. Demonstrating 
that such naive procedure leads to a problem, when the continuum limit is 
taken, we will introduce an improved scaling test. By also discussing an al- 
ternative scaling test, we will then demonstrate that physical results and the 
approach to the continuum limit can strongly depend on the choice of the 
setup for a scaling test. We provide an explanation for this finding by the 
effects of exact chiral zero modes of the overlap Dirac operator taken in the 
valence sector. 



4-.1 Naive Scaling Test 



The naive scaling test follows very closely the example of free fermions dis- 
cussed above. As in the free case, we match the pseudoscalar meson masses 
computed from the correlator constructed from the pseudoscalar interpolat- 
ing field (PP correlator, Cpp(t)), which appears to be the most natural choice. 
Of course, in practical simulations, the pseudoscalar masses cannot not be 
matched exactly. We therefore use that value of the overlap pseudoscalar mass 
which is closest to the targeted twisted mass value. The corresponding over- 
lap bare quark mass then defines our matching quark mass. Using Cpp(t) for 
the matching condition we will refer to as naive matching condition. To pro- 
ceed, we evaluate the pseudoscalar correlation function Cpp(t) for light quark 
masses using ensembles -B^^vi, Q and Di. 



The matching procedure is illustrated for the case Bey^ in Fig. 3(a) The 



values of the matching quark and pseudoscalar masses obtained from different 
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Ensemble X 


arrin 




/7 'yy~t m^tch 


"'"'ov 


matching — > 




naive 


improved 


alternative 


Bh 


0.1961(17) 


0.015 


0.016 


0.017 


Be,Vi 


0.1592(24) 


0.007 


0.011 


0.009 


Bl.V2 


0.1389(14) 


0.007 


0.009 


— 


Be,V3 


0.1359(7) 


0.008 


0.008 


— 


Ch 


0.1520(15) 


0.011 


0.012 


0.013 


Ci 


0.1209(40) 


0.005 


0.006 


0.006 


Dh 


0.1252(13) 


0.009 


0.010 


0.010 




0.0980(19) 


0.002 


0.004 


0.004 



Table 2 

The pseudoscalar masses and overlap matching quark masses for different ensembles. 
Results of employing different matching conditions are shown. 



ensembles are tabulated in table [2l In Fig. 3(b) we show the pseudoscalar 
decay constant computed from the valence overlap fermions at the matching 
mass as a function of a^. As in the case of free fermions, we indeed find 
a scaling of the pseudoscalar decay constant compatible with the expected 
0{a^) behaviour. However, we encounter a problem with this naive scaling 
test when extrapolating the mixed action pseudoscalar decay constants to the 
continuum limit, we find ro/°^ = 0.236(9) while the continuum extrapolated 
value for unitary twisted mass fermions comes out to be tq/*™ = 0.181(10). 



Thus, it seems that we are obtaining an inconsistent continuum limit from our 
mixed action simulation. Of course, this observation demands an explanation, 
since this is clearly in contradiction to our expectation and must be an artifact 
of our setup. Indeed, the setup of our naive scaling test is too naive. It does 
not take into account the above mentioned mismatch of the eigenvalue spectra 
of the lattice Dirac operators. In particular, the pseudoscalar correlator might 
be strongly affected by zero modes of the overlap Dirac operator, especially in 
the small volume and at the small value of the quark mass employed by using 
ensembles -B£,Vi, Ci and De. 

The suspicion of the overlap zero modes being responsible for the apparent 
inconsistent continuum limit values for the pseudoscalar decay constant leads 
to an improved scaling test. 
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Fig. 3. (a) Naive matching condition for the case of ensemble -Bf^Vi- (b) Continuum 
hmit of /,r in the case of the naive matching condition. 



4-2 Improved Scaling Test 



In order to check the assumption that the zero modes play a special role in 
our mixed action setup, we should use observables for which such zero mode 
contributions are cancelled and which are therefore in this sense "improved" . 
Such observables can then be used both in the matching condition and to com- 
pute the physical quantity. As we will demonstrate below, once an improved 
version of -in our case- /,r is used, all matching conditions considered lead to 
completely consistent continuum limit values. 

Our improved scaling test consists of studying instead of the pure pseudoscalar 
correlation function the correlator 



Cpp_ss(t) = Cpp(t)-Css(t), (6) 

where the scalar correlation function is subtracted and hence the contribution 
of zero modes is exactly cancelled, since they are eigenf unctions of 75. In this 
way, we match the overlap (mixed action) pseudoscalar meson mass resulting 
from such correlator to a unitary pseudoscalar meson mass extracted from 
the standard PP correlator. This defines the improved matching condition. 
Moreover, we extract from the PP-SS correlator. Therefore, we have an 
improved setup in both the matching condition and in the observable itself - 
in both we avoid the effects of zero modes. 



We show the matching of the pseudoscalar mass in Fig. 4(a) and the scaling in 



of the pseudoscalar decay constant in Fig. 4(b) where we again observe the 
expected O(a^) scaling behaviour. Performing a continuum limit extrapolation 
of /tt, we find ro/°^ = 0.196(13) which is now completely compatible with the 
one of twisted mass fermions, Tq/*™ = 0.181(10). 
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.04 



(a) (b) 

Fig. 4. (a) Improved matching condition for the case of ensemble B^y-^. (b) Contin- 
uum hmit of /tt in the case of the improved matching condition. 




Fig. 5. (a) Matching of overlap and MTM quark masses (ensemble -B£,Vi) at a heavy 
quark mass value (yielding rQiri-,^ = 1.5, which is denoted by vertical line), where the 
zero modes effects are strongly suppressed. Idea of the matching is explained in the 
text, (b) Continuum limit of f-,^ in the case of the alternative matching condition. 

4-3 Alternative Scaling Test 



The results of the previous two sections indicate that indeed the zero modes 
play a special role in the mixed action setup used here. Using the correla- 
tor Cpp_ss(^) which cancels the zero mode contributions exactly is not the 
only possibility to reach this goal. Our alternative scaling test starts by the 
observation that zero modes can also be suppressed by using either a large 
enough volume or large quark masses, see also the discussion below. One way 
of avoiding the effects of chiral zero modes, at least in the matching condition, 
is then to match the theories at a heavy quark mass values where the effect 
of the zero modes is strongly suppressed. 

An idea of such matching is the following. We choose a value of rom^ = 1.5 
in the regime of quark masses where the difference between the C-p-p{t) and 
Cpp-ss(^) correlators is negligible, which indicates that the effects of the zero 
modes are small. The dependence of the pseudoscalar mass on the overlap 
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Fig. 6. Continuum limit scaling of ?"o/7r, using different matching conditions and 
extracting from the PP-SS correlator. All matching conditions lead to a consistent 
continuum limit with the one of the unitary approach. In the graph the data for 
"naive matching" is shifted to the left and the data for "alternative matching" to 
the right, for reasons of a better presentation. 




and MTM quark mass (using partially quenched MTM data) allows us to find 
the values of quark masses which correspond to the chosen value of rom.^. 



As can be seen on Fig. 5(a), the smallness of the zero modes effects at this 
value of romir means that the difference /'o(/°^ — /^™) is significantly reduced 
with respect to the large difference in ro/,r in the small quark mass regime. 
The ratio of MTM and overlap bare quark masses that lead to rom,r = 1.5 
provides an estimate of the ratio Rzp = Zp^/Zp™. In this way, we can define 
the alternative matching mass as RzpfJ^q, where fig is the unitary MTM quark 
mass. The values that we find from this procedure (for all ensembles with 
linear lattice extent L ^ 1.3 fm) are again tabulated in table [2j 



Once we have the matching masses from the alternative matching condition, 
we can perform the continuum limit scaling of the corresponding decay con- 
stant, extracted again from the Cpp_ss(^) correlator, such that the zero modes 
effects are suppressed both in the matching condition and in the values of /yr 
(in this sense, such setup is again an improved one with respect to the naive 
scaling test). Needless to say that we again observe a nice O(a^) scaling of the 
pseudoscalar decay constant towards the continuum limit (Fig. 5(b)), where 
we find ro/°^ = 0.203(19), which is to be compared to the continuum limit 
value of twisted mass fermions ro/*™ = 0.181(10). Again, a consistent result 
is found. Note that a similar conclusion is obtained if instead of performing 
the alternative matching with partially quenched data, we use the unitary 
ensembles Bh, Ch and Dh, at the heavier masses, to impose the matching 
condition. 
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5 Interpretation and Discussion 



In the previous section, we have discussed continuum limit scaling tests, em- 
ploying three different matching conditions and using two definitions for /°^. 
The disturbing observation in these scaling tests has been that in case of the 
naive scaling test, an apparent inconsistency of the continuum limit value of 
has been encountered which, however, could be resolved by employing an 
improved and an alternative scaling strategy. One may wonder then, whether 
in the naive scaling test the matching condition, the used observable or even 
some interplay between both are responsible for the inconsistent continuum 
limit of /,r- 

Fig. E] provides an answer to this question, at least in the case of considered 
here. In the graph, we use the naive, improved and alternative matching con- 
ditions but evaluate f°^, in all cases, from Cpp_ss(^)- Obviously, with all three 
matching conditions we obtain a completely consistent continuum limit value 
of /°^, agreeing well with the one from the unitary setup. Thus, it seems that 
the choice of the observable has a significant impact on the observed mismatch 
in the naive scaling test of /,r due to the zero modes effects. 

It is interesting to also perform the opposite test by using again all three 
matching conditions but this time only use computed from the PP cor- 
relator. In this case, we obtain inconsistent continuum limit values, ro/°^ = 
0.236(9), ro/;^ = 0.252(10) and ro/^^ = 0.256(10) for the naive, improved and 
alternative matching conditions, respectively, while for the unitary setup we 
have the value tq/*™ = 0.181(10). As a consequence, it is important to define 
observables such that they are not influenced by zero modes. 

As we have shown above, mixed actions involving chiral fermions only in the 
valence sector can exhibit particular effects. Indeed, the spectra of overlap 
and twisted mass Dirac operators differ in an essential way. The overlap Dirac 
operator, being chirally symmetric, incorporates the physical effects of exact 
chiral zero modes at any value of the lattice spacing. This is not the case 
for the Wilson twisted mass Dirac operator. In the context of a mixed action 
this implies that the effect of zero modes in the valence sector is not properly 
suppressed by a fermionic determinant which depends only on sea quarks. Zero 
mode effects can therefore introduce potential difficulties when comparing the 
continuum limit extrapolation of mixed action data with respect to the unitary 
case. 

The effects of zero modes are observable dependent. In particular, when aiming 
at studies of QCD on a large volume, an appropriate choice of the interpolating 
operator can be used to reduce or to cancel the effects of zero modes on a 
given observable. Indeed, zero modes effects can be interpreted as finite size 
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effects and are thus suppressed as the volume is increased. Moreover, when 
the quark mass increases, the relative contribution of the zero modes decreases 
with respect to the one of the other low-lying modes. Observables receiving 
contributions from zero modes are therefore less influenced by their effects as 
the quark mass increases. 

Before closing this section, we collect a few remarks about the continuum limit 
scaling of the mixed and unitary actions. 

• The previous sections and discussion have illustrated the influence of ex- 
act zero modes of the overlap operator in a mixed action setup. We have 
examined three examples of continuum limit scalings which give partially 
different results. Of course, our list is by no means exhaustive, nor do we 
claim to have provided the "optimal" continuum limit strategy. However, 
we have clearly provided a warning that in a mixed action setup special 
care has to be taken when lattice Dirac operators with exact chiral zero 
modes are employed in the valence sector. In the following, we will make 
an attempt to provide a quantitative estimate in the region of volume and 
pseudoscalar mass where the zero modes are the dominant effect and where 
they can be safely neglected. 

• The Cpp_ss(t) correlator in eq. (|6]) has different discretization effects than 
Cpp(t). Indeed, it is known that the scalar correlator Css(t) is especially 
vulnerable to the double pole contribution to the meson propagator [22]. 
This unphysical effect is present in all unitarity violating setups (quenched, 
partially quenched and mixed action) and a separate analysis indicates that 
the scalar correlator indeed obtains a negative contribution of this kind, 
which has opposite sign with respect to the one from zero modes |19]. 

• The Cpp-ss(^) correlator receives excited state contributions from the scalar 
correlator. This can affect the extraction of the pseudoscalar ground state 
in the large quark mass region. In the regime of quark masses that we 
consider in this work this appears to be a negligible effect. In general, we 
observe clear signals for the plateaux, which we have furthermore controlled 
by comparing the extracted values of m^r and to an analysis including 
an estimate of the systematic effect due to the choice of the fit interval, see 

also ng. 

• For twisted mass fermions, the use of Cpp_ss(t), instead of Cpp(t), has a 
negligible impact on the values of /^r and m-Tr in the regime of parameters 
(quark mass, lattice size and lattice spacing) that we consider in this work. 
We only consider the latter case. 

• In the Wilson twisted mass formulation, a potential contribution from the 
large discretization effects present in the neutral pion mass can occur, in 
particular through finite volume effects. However, these cut-off effects are in 
principle properly taken into account by the fact that both the quark mass 
and the volume are kept fixed during the continuum limit extrapolation. 
The fact that we observe the expected C(a^) scaling behaviour, gives us 
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Ensemble 


unitary MTM 


rn t-rr 
' U J TT 

naive matching 


rn f-rr 
' U J TT 

improved matching 

X X kJO V.^"^! 1 ClClUJi 


Bh 


3413f33') 


3973f65l 


3626fl07) 


Ch 


0.3373(33) 


0.3774(63) 


0.3669(78) 


Dh 


0.3382(50) 


0.3493(100) 


0.3558(90) 


continuum limit 


0.334(11) 


0.330(17) 


0.345(18) 



Table 3 

The values of the pseudoscalar decay constant for ensembles -B/j, and at the 
matching mass (with the naive and improved matching condition). The continuum 
limit values are also shown. 

confidence that potential isospin breaking effects are indeed removed in the 
continuum limit, as also observed in [36] . 



6 Quark Mass Dependence and Finite Size Effects 

In this section, we will perform a matching of the twisted mass sea and valence 
overlap actions at a heavier quark mass and at larger volumes. The purpose 
of this investigation is twofold. First of all, these investigations will tell us on 
a quantitative level at which value of the quark mass and volume the effects 
of the exact chiral zero modes are harmless (safe region), when their effect 
is significant (hazardous region) and when their effect is unacceptably large 
(non-safe region). Second, using larger quark masses and volumes serves as 
a test of the picture developed in the previous sections about the role of the 
zero modes in a mixed action simulation. 

6.1 Quark Mass Dependence 

For the matching of the heavier quark mass we use the ensembles Bh, Ch and 
Dh of tabled! We apply for these ensembles both the naive matching procedure 
(naive matching condition using from the PP correlator) and the improved 
one (improved matching condition using from the PP-SS correlator). In 
both cases we observe, as expected, a nice scaling of the pseudoscalar decay 
constant in towards the continuum limit. In table |3] we give the values of 
the pseudoscalar decay constants from the naive and the improved matching 
conditions, as well as one for the unitary twisted mass case. As can be seen in 
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the table, the differences in the values of and /*™ for all ensembles are much 
smaller than in the case of Biy^ , Ci and Di ensembles and all continuum limit 
values are consistent with one another. This clearly confirms our interpretation 
that the zero modes effects are large at a quark mass corresponding to a 
pseudoscalar mass of 300 MeV (in infinite volume). At 450 MeV, these effects 
are strongly suppressed although not totally negligible (see table |3]). 



6.2 Finite Size Effects 



Finite size effects (FSE) in the light pseudoscalar meson observables consid- 
ered in this work can have multiple origins. As mentioned earlier, zero modes 
of the Dirac operator can introduce significant finite volume effects in the pseu- 
doscalar correlator. Other sources of FSE are related to pion loops or to the 
actual size of the box with respect to the typical hadronic scale. Concerning 
pion loops, in the case of twisted mass fermions a contribution to FSE from the 
lighter neutral pion mass can in principle be present at finite lattice spacing, 
although this is expected to vanish when performing the continuum limit ex- 
trapolation (as discussed in the previous section). When considering the scalar 
correlator, an additional FSE arises from the bubble diagram describing the 
exchange of two pseudoscalar mesons. Isolating the relative contribution of 
these effects in our data for the mixed and unitary actions is beyond the scope 
of this work. Here, we rather concentrate on the comparison of the dependence 
of and on the volume, since this should give some evidences of the 
presence of FSE from the topological zero modes. 

We consider three ensembles Biy^^, Biy^, Biy^ at the coarsest value of the 
lattice spacing, see tabled! The corresponding linear lattice extents are 1.3 fm, 
1.7 fm and 2.0 fm, respectively. The quark mass is set to rriq ~ 20 MeV and 
the values of m^rL at the matching point are 2.5, 2.8, 3.3, for ensembles B^y^, 
Biy^, Biy^, respectively. 

To illustrate the magnitude of finite size effects, we plot the quark mass depen- 
dence of the pseudoscalar meson decay constant - in this section we only use 



the PP correlator - for our three ensembles in Fig. 7(a) The horizontal bands 
correspond to the unitary values of /*™. The matching mass (employing the 
naive matching procedure) is shown by the vertical dashed line. We observe 
significant FSE for both the mixed and unitary actions when reducing the 
lattice size from 1.7 fm down to 1.3 fm, while they are smaller between 2.0 fm 
and 1.7 fm. Note that for there is a relative sign in the contribution to FSE 
coming from pion loops and zero modes. This could induce, in practise, some 
cancellations of FSE in the mixed action data. 

The differences between the decay constants and /*™, at the matching 



18 



L/a=1 6 Overlap > 
MTM ■ 

L/a=20 Overlap 
MTM 

L/a=24 Overlap 
MTM ■ 



0.005 0.01 0.01 5 0.02 0.025 0.03 0.035 



0.001 0.002 0.003 0.004 0.005 



(a) 



(b) 



Fig. 7. Finite size effects from ensembles Biyi, Bpy^, B^y^ (see tabled]), (a) Quark 
mass dependence of the pseudoscalar meson decay constant, (b) Volume dependence 
of the difference between the pseudoscalar meson decay constant in the mixed and 
unitary actions. We observe that the difference between /*™ and decreases when 
increasing the volume. 



mass, are plotted in Fig. 7(b) In such difference, many of the contributions 
to FSE are expected to cancel. The fact that — /*™ has a non vanishing 
slope when changing the lattice size can be interpreted as a FSE coming from 
the zero modes of the overlap operator. This hypothesis is supported by the 
fact that this slope is within errors compatible with zero when instead of 
using /°^[PP] in the difference, one considers /°^[PP — SS], extracted from 
the C-p-pssif) correlator. 



Fig. 7(b) provides indications that FSE from topological zero modes are not 
totally negligible for the ensemble Bg^y^, with our largest lattice size L ^ 2.0 fm 
and tUt^L ^ 3.3. For the heavier ensemble B^, with L ^ 1.3 fm and m^L ^ 3.1 
and the same value of the lattice spacing, we observed that af°^ — af^ = 
0.011(1), arising mainly from unitarity violations, since zero mode effects were 
observed to be small in this case. By taking this value as a rough estimate for 
the size of unitarity violations at this value of the lattice spacing, we estimate 
that a lattice size L > 2.4 fm (corresponding to itLt^L > 4) would be needed, 
at the lighter quark mass rrig ~ 20 MeV, to enter in the regime of small FSE 
from chiral zero modes. In practise, this would require to deal with the overlap 
operator on an L/a = 32 lattice. 

We close this section by a remark about the size of the unitarity violations in 
our mixed action. The difference — /*™ in Fig. 7(b) is numerically large 



for all three volumes. This difference is not expected to vanish in the infinite 
volume limit due to the presence of unitarity violations in the mixed action 
data at non-zero values of the lattice spacing. At a value of the lattice spacing, 
a ^ 0.08 fm, with L ^ 2.0 fm and m^rL ^ 3.3, we observe that the mixed and 
unitary values of fj^ can differ by 20%. Fig. [6] illustrates a similar effect on a 
smaller volume. Indeed a significant difference between overlap and twisted- 
mass data is observed at the coarser value of the lattice spacing, a ~ 0.08 fm. 
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This points towards the necessity of a proper control of the continuum hmit 
extrapolation when working in a mixed action setup. 



7 Conclusion 

We have studied the continuum limit scaling of a mixed action with overlap 
valence quarks on a twisted mass sea employing three values of the lattice 
spacing. We observe the expected O(a^) scaling behaviour for the here in- 
vestigated pseudoscalar decay constant. We have demonstrated that in such 
mixed action approach the zero modes of the overlap operator which are un- 
matched by the twisted mass operator at the lattice spacings employed here, 
can play a special role. They can severely affect physical quantities such as 
the pseudoscalar decay constant. 

The precise values of the quark mass and/or lattice size at which the effects 
of zero modes become relevant for a given observable are not known a priori. 
A dedicated numerical study is therefore needed to address this question. To 
this end, we performed simulations at two values of the quark mass and three 
volumes. These calculations then provided information on the values of the 
parameters for which these zero mode effects are significant or when they 
can be neglected. Our findings are summarised in Fig. [HI There we indicate 
on a more qualitative level the regions in the quark mass and volume where 
simulations ought to be safe, i.e. where the zero modes play no significant 
role, the hazardous regions, where their effect can be significant and a non- 
safe region, where their effect become really large. 

We want to stress that the difficulties in the continuum limit approach of /^r 
arising from the zero modes of the overlap Dirac operator originate solely from 
the particular definition of /^r and not from the matching condition employed. 
We demonstrated this clearly by, employing three different matching condi- 
tions, but using only a definition of (constructed from Cpp_ss(^)) which 
is not affected by zero modes, see Fig. O On the contrary, by employing a 
definition of which is highly sensitive to zero modes (i.e. constructed from 
Cpp(t)), we find inconsistent continuum limit results for all three matching 
conditions employed. This might suggest that if it is possible to construct 
physical observables such that they are not affected by zero modes consis- 
tent continuum limit values can be obtained. However, this expectation needs 
further corroboration by studying other physical quantities. 

The effects of zero modes depend on the choice of the operator. Furthermore, 
their effects vary from one observable to another. We refer to a forthcoming 
publication |19] for a study of the impact of zero modes on other observ- 
ables and operators. For some observables, isolating the contribution from 
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Fig. 8. A qualitative graph indicating the safe, hazardous and non-safe regions 
of parameters in physical units (linear extent of the lattice vs. pseudoscalar meson 
mass) in mixed action simulations with overlap valence and twisted mass sea quarks. 
The parameter values in the "safe" region are such that the effects of chiral zero 
modes of the overlap operator are negligible. The boundary between the "safe" and 
"hazardous" region corresponds to m-,^L = 4 and the boundary between "hazardous" 
and "non-safe" is m-j^L = 3. 

zero modes can be a difficult task. For the study of such observables, it is 
particularly important to perform the simulations in a regime of parameters 
where these effects are negligible. 

By working in the safe regime, where the effects of zero modes are negligible, 
we believe that this mixed action approach should allow to obtain accurate de- 
terminations of quantities for which good chiral properties of valence fermions 
are essential, such as the kaon bag parameter or the amplitudes of — t- tttt 
decays. 
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